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ABSTRACT 

Within the frame of the NUclei of GAIaxies (NUGA) project, we have determined the distribution and kinematics of the molecular gas within 
^ the central kpc with high spatial resolution (100-150pc), for a sample of active galaxies. The goal is to study the gas-fueling mechanisms in 
l AGN. 

CO t We present interferometric observations of '-CO(l-O) and '^CO(2-l) line emission from the Seyfert 2 galaxy NGC 6574, obtained with the 
IRAM Plateau de Bure Interferometer (PdBI). These data have been combined with 30m mapping data in these lines to correct for the flux 
resolved by the interferometer. At an angular resolution of 0.7" (= 1 lOpc), the '^C0(2-l) emission is resolved into an inner disk with a radius 

C\i of 300 pc. 

T-H The molecular gas in NGC 6574 is primarily distributed in four components: nucleus, bar, spiral arms - winding up into a pseudo-ring - and 
C an extended underlying disk component. For the overall galaxy host, we find a '^CO(2-l) to '^CO(l-O) line ratio of ~0.3 indicative of cold or 
sub-thermally excited gas. For the nucleus, this ratio is close to unity, indicating emission from dense and warm molecular gas. Modeling the 
^ gas kinematics with elliptical orbits shows that the molecular gas in the differentially rotating disk of NGC 6574 is strongly influenced by the 
• presence of a stellar bar. The nuclear component shows an extension toward the southeast that may be an indication of the lopsidedness of the 
^\ nuclear gas distribution. 

We computed the gravity torques exerted from the stellar bar on the gas, deriving the gravitational potential from near-infrared images, and 
weighting the torques by the CO distribution. We find negative torques for the gas inside the ring, since the gas aligned with the bar has a slight 
advance phase shift, leading the bar. This means that gas is flowing in towards the center, at least down to 400pc in radius, which can explain 
the observed high nuclear gas concentration. This concentration corresponds to a possible inner Lindblad resonance of the bar, according to the 
measured rotation curve. The gas has been piling up in this location quite recently, since no startburst is been observed yet. 

Key words, galaxies: individual: NGC 6574 - galaxies: active - galaxies: kinematics and dynamics 



1. Introduction 

The study of interstellar gas in the nuclei of galaxies is a 
fundamental for understanding nuclear activity through fu- 
eling the active galactic nuclei (AGN) and their relation to 
circum-nuclear star formation. Some studies claim that pure 
gaseous density waves (spirals, bars, warps, or lopsided insta- 
bilities) may be driving gas inflow towards the AGN (Heller 
& Shlosman 1994; Elmegreen et al. 1998; Regan & Mulchaey 
1999). While such dynamical perturbations are responsible for 
the infall of gas on large scales, the processes responsible for 



Send offprint requests to: A. Eckart; eckart@phl.uni-koeln.de 



removing angular momentum on small scales (sub-kpc) have 
not been understood very well yet. 

Different scenarios have been introduced trying to explain 
this phenomenon, such as nested bars (e.g., Friedli & Martinet 
1993, spirals e.g., Martini & Pogge 1999,) warped nuclear 
disks (Schinnerer et al. 2000a,b), and lopsidedness or m=l 
instabilities (Kormendy & Bender 1999; Garcia-Burillo et al. 
2000). Most of the gas within the central kiloparsec of spiral 
galaxies is in the molecular phase, while the atomic gas is de- 
ficient there. 

This makes CO lines the best tracers of nuclear gas dynam- 
ics in the nuclear interstellar medium. High-resolution (l"-2") 
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maps of the molecular component in the centers of galaxies are 
required to model these nuclei within the NUclei of GAlaxies 
(NUGA) project. The NUGA project aims at establishing a 
high angular resolution 0.5" - 1") and high-sensitivity CO 
survey of twelve objects. It covers the whole sequence of activ- 
ity types: Seyfert 1, Seyfert 2, LINERs, and transition objects, 
and it comprises 12 objects as a core sample and a total of about 
30 objects as an extended sample. The survey is being carried 
out with the IRAM Plateau de Bure mm-Interferometer (PdBI) 
in France, because it offers the best combination of sensitivity 
and resolution, both crucial for this project. A more detailed 
description of the NUGA project is given in Paper I by Garcia- 
Burillo et al. (2003). 

Previous surveys of molecular gas in active galaxies have 
been carried out by other groups (Heckman et al. 1989; 
Meixner et al. 1990; Vila-Vilaro et al. 1998; Sakamoto et al. 
1999a,b), but have had insufficient spatial resolution (4-7") 
to resolve the nuclear disk structures, or else were limited to 
small samples (Tacconi et al. 1997; Baker 2003). Besides case- 
by-case analyses and simulations of each object of the sam- 
ple (Garci'a-Burillo et al. 2003; Combes et al. 2004; Kiips 
et al. 2005; Garci'a-Burillo et al. 2005), the collected data of 
the NUGA project also provides an initial statistical basis for 
studying different mechanisms that may be responsible for gas 
flow toward the nuclei or may account for further accretion in- 
ward. 

This paper describes the distribution and dynamics of 
molecular gas in NGC 6574, one of the galaxies belonging to 
the NUGA survey. NGC 6574 is a Seyfert 2 galaxy of Hubble 
type SB(s)bc, at a distance of 33 Mpc. This results in a spa- 
tial scale of ~160 pc per arcsecond. Its main characteristic is 
a symmetric structure with two spiral arms to the north and 
south of the center, nested on a small bar within the central 8". 
Kotilainen et al. (2000) present high spatial-resolution, near- 
infrared broad-band JHK and Bry 2.166 /im and H2 1-0 S(l) 
2.121 yum emission line images of the circum nuclear star- 
formation ring. The overall near-infrared and radio morpholo- 
gies (see Laine et al. 2006) generally agree with each other. The 
observed H2/Bry ratio indicates that the main excitation mech- 
anism of the molecular gas is UV radiation from hot young 
stars. Shocks are likely to contribute only in a few regions. To 
explain the NIR data, Kotilainen et al. (2000) prefer the model 
of an instantaneous burst of star formation with an upper mass 
cutoff of M„ - 100 M0 occurring ~6-7 Myr ago. An analy- 
sis of the molecular gas distribution at an angular resolution of 
about 4" in the '^CO(l-O) line is presented by Sakamoto et al. 
(1999ab). They find a prominent CO peak of 1 kpc diameter at 
the center of the galaxy on a smooth gas disk with two gas/dust 
lanes that extend about 2 kpc from the central peak to east and 
west at a PA of ~105°. From the tips of this gaseous bar, spiral 
arms start to form a ring-like structure. 

In this paper we analyze data obtained with the IRAM 
PdBI and the 30 m single-dish telescope in the '^CO(l-O) and 
'^CO(2-l) lines. The angular resolution we obtained in the 
combined data set is about 2" and 0.7". We describe the molec- 
ular gas distribution and kinematics with elliptical bar orbits, 
to account for the influence of the stellar bar potential in the 
nucleus of NGC 6574. In Section 2 we describe the CO obser- 



vations and in Section 3 we present the first results. The kine- 
matic model and its results are discussed in Section 4. Gravity 
torques and the nature of gas flows are derived in Section 5, 
while Section 6 concludes. 



2. Observations and data reduction 

In this section we present the '^CO observations with the 
IRAM 30m single dish, with the IRAM PdBI, including the 
short-spacing correction. 



2.1. ^^CO observations with the IRAM 30m single dish 

We performed IRAM 30m observations in a 7 x 7 raster pat- 
tern with 10" spacing in August 2005, in order to add the short 
spacings and to estimate the flux density filtered out by the in- 
terferometric observations. We used 4 SIS receivers to observe 
simultaneously at the redshifted frequencies of the '^CO(l-O) 
and the '^CO(2-l) lines, thereby assuming the LSR velocity 
of NGC6574 of 2282 km s '. For each fine, a 4 MHz chan- 
nel filter bank was used with a velocity resolution smoothed 
to about 10 km s"' in order to improve the signal-to-noise 
ratios of the individual spectra. The 4 MHz filter bank con- 
sists of nine units. Each unit has 256 channels which covers 
a total bandwidth of 1 GHz. At 115 GHz and 230 GHz, the 
HPBW is 21.5" and 10.75", respectively. The parameters for 
the two observing runs are summarized in Table [T] We used 
the wobbler switch mode with a 6 s cycle and a throw of 240". 
With the forward efficiency Fe/f, the main-beam efficiency is 
Beff^FeffXrjT^^ = 0.74 at 115 GHz and 0.53 at 230 GHz. 
The typical system temperature varied between 200 and 450 K 
at both frequencies. The pointing was checked regularly on 
continuum sources and the accuracy was consistent to 3" rms 
achieved by Greve et al. (1996). The total bandwidth available 
was 512 MHz at 1 15 GHz and 230 GHz, corresponding to 1330 
km s ' and 666 km s respectively. The data reduction and 
analysis were done with the GILDAS package (e.g. Guilloteau 
& Lucas 2000). For the short-spacing correction and further 
investigations, only the 4 MHz filter bank was used. 



Parameters 


CO(l-O) 


CO(2-l) 


central coordinates 


RA = 18'' 


ir"51.3', DEC = -i-14°58'52.l" 


backends 


4 MHz filter bank/autocorrelator 


redshift z 




0.007662 


sky frequency [MHz] 


1 14.400 


228.798 


beam size [arcsec] 


21.50 


10.75 


beam efficiency B^ff 


0.739 


0.527 


forward efficiency F^ff 


0.95 


0.91 


total obs. time [sec] 


14994 


14994 


Nr. of positions 


49 


49 


Nr. of scans 


294 


294 



Table 1. Parameters of the observing run with IRAM 30m sin- 
gle dish. 
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2.2. ^2 CO observations with IRAM PdBI 



2.3. Short-spacing correction 



We observed the emission of the J=l-0 and J=2-l lines of 
'^CO in a single field centered at the radio position of the 
AGN (i.e., Q'2000=18*'ll™ 51.3^and 52ooo=14° 58'52.1") in 
NGC 6574 with the IRAM Plateau de Bure interferometer 
(PdBI) in April 2002 (D configuration), in November 2003 (C 
configuration) and February 2004 (A and C configurations). 
The six 15 m antennae were equipped with dual-band SIS re- 
ceivers. The spectral correlators were centered at 1 14.400 GHz 
and 228.798 GHz, respectively, with three correlator units cov- 
ering a total bandwidth of 400 MHz at each frequency. The dif- 
ference between LSR and heliocentric velocities is 12 km s 
therefore the observations were centered on Vhei - 2282 km s"' . 
The correlator was regularly calibrated by a noise source in- 
serted in the IF system. 

The bandpass calibrator was 3C 273 or MWC349, while 
the phase and amplitude calibration were performed on the 
nearby quasars 1923+210 and 1749+096, depending on the 
data set. Flux densities were calibrated relative to CRL618 and 
MWC349. The data were phase calibrated in the antenna-based 
mode. The frequencies were centered on the redshifted '^CO(l- 
0) line in the USB at 3 mm and on the redshifted '^€0(2-1) 
line in the LSB at 1 mm. For each line, the total bandwidth was 
400 MHz and the spectral resolution 1 .25 MHz. The integration 
time for the central pointing amounts to ~16 hrs on the source. 
The water vapor ranged between 4 and 10 mm (i.e., opacities 
of ~0.2-0.3) resulting in system temperatures of approximately 
200-300 K on average. 

The flux densities of the primary calibrators were deter- 
mined from IRAM measurements and taken as input for de- 
riving the absolute flux density scales for our visibilities, esti- 
mated to be accurate to 10%. The parameters for the four ob- 
serving runs are summarized in Table |2] 

The data reduction and mapping was performed using the 
GILDAS software. Data cubes with 512 x 512 spatial pixels 
(0.25"/pixel) were created with velocity planes separated by 

5 km s ' (for PdBI alone, and by 10 km s ' for the short- 
spacing corrected data). The cubes were cleaned with the Clark 
(1980) method for ^^CO{l-Q) data and with MX for '^00(2-1) 
data. As synthesized clean beams we used a 2.48 " x 1.69" 
Gaussian with PA = 26° at 3 mm and a 0.93 " x 0.5" Gaussian 
with PA = 12° at 1 mm. The rms noise levels in the cleaned 
maps (at 5 km s"' velocity resolution) are 3 mJy/beam and 

6 mJy/beam for the '^00(1-0) and '^00(2-1) Unes, respec- 
tively. No continuum emission was detected toward NGC 6574. 
The maps were corrected for primary beam attenuation. 



Parameters 



central coordinates 
redshift z 

sky frequency [MHz] 
beam size [arcsec] 
position angle [degree] 



CO(l-O) 
RA= IS"!!" 

114.400 
2.48x1.69 
26 



CO(2-l) 



51.3', DEC = +14°58 52.1 
0.007662 

228.798 
0.91x0.48 
12 



Short spacings were included using the SHORT-SPACE task in 
the GILDAS software. To find the best compromise between 
good angular resolution and complete restoration of the miss- 
ing extended flux, the weights attached to the 30m and PdBI 
data were varied. At the end, a factor of two was applied to 




00 200 



velocity in kiWs 



CO(2-l) 




Table 2. Parameters of the observing run with IRAM PdBI. 



velocity in km/s 

Fig.l. The '^CO(l-O) and '^00(2-1) line intensity derived 
from the PdBI data set (red) and the combined PdBI+30m 
data set (black). The spectra combine the information shown 
in the channel map in Fig.|2] The comparison shows the larger 
amount of extended emission for the '^CO(l-O) line that is re- 
solved by the interferometer. Both sets of spectra also demon- 
strate the asymmetry in flux density between the red and blue 
shifted part of the spectrum. The red side of the spectrum is 
brighter mostly in its extended emission. 



the weights of the 30m data enabling us to recover more than 
40% of the missing flux and to retain the angular resolution of 
2.48" X 1.69" at PA of 26° for ^^CO{\-Q) emission and 0.93" 
X 0.5" at PA of 12° for '^00(2-1) emission. The combined 
data sets were written to visibility tables, converted to maps 
using standard reduction procedures, and then deconvolved us- 
ing the Clark algorithm. The weights were adjusted to get the 
same mean weights in the single-dish data as in the interferom- 
eter data in the u-v range of 1.25 DjA to 2.5 D//1(D=15 m). 
All figures presented in this paper are done with short-spacing 
corrected data. 

The line intensities and CO luminosities integrated over the 
entire linewidth of the full Av=370 km s ' extent of the com- 
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R.ii 

Fig. 2. '^CO(l-O) velocity-channel maps observed with the PdBI and corrected through short-spacing by the 30m single dish in 
the nucleus of NGC 6574. The clean beam resulting from the observations is 2.48" xl.69" at PA=26°. It is plotted as a filled 
ellipse in the bottom left corner of each panel. We show a field of view of 55". The phase-tracking center is indicated by a cross 
at ©2000 = 18''11'"51.3' and 620m = +14°58'52.1". Velocity channels are displayed from v = -180 km s^' to v = 230 km s^' in 
steps of 10 km s"^ Velocities are in LSR scale. Contour levels are 0.42cr, 0.84cr to 4.8cr in steps of 0.42cr, where the 1-cr rms is 
cr=2.39 mJy beam"'. 



bined PdBlH-30m map after correcting for primary beam re- 
sponse are listed in Table[3] Since these maps contain the entire 
galaxy, the corresponding quantities give information on the 
overall molecular excitation in NGC 6574. The difference be- 
tween the integrated quantities for the '^CO(l-O) line emission 
already indicate the presence of an extended disk component 
that is resolved out by the interferometer (see Fig.[T]). 



3. Results 

In this section we first describe the morphology of NGC 6574 
and show its kinematics, then estimate the gas masses and re- 
lated parameters. 

3.1. Morphology of the CO emission 

In Fig.|2]we show the velocity-channel maps of the ^^COCl-O) 
line emission in the central region of NGC 6574, with a veloc- 
ity range of 410 km s"' and a velocity resolution of 10 km s^'. 
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relative R,A, in arcsec 

Fig. 3. Primary beam corrected integrated intensity maps for the '^CO(l-O) (left) and '^CO(2-l) (right) line emission. The straight 
lines at PA - 165° indicate the position of the major axis of NGC 6574 along and perpendicular to which we obtain position- 
velocity diagrams. The '^CO(2-l) map is shown in the same frame as the '^CO(l-O) map. The dashed circles outline the central 
regions that approximately include the FWHM areas of the PdBI primary beams at the two frequencies. For larger radii the map 
is too noisy. Inside that radius the nuclear emission as well as emission from the bar can be detected. This includes the tips of 
the bar at which the spiral arms start. Contour lines are 1, 2, 3, 4, .... 16 Jy km/s /beam for '^CO(l-O) and 1.8, 3.6, 5.4, 7.2, ... 15 
Jy km/s /beam for '^CO(2-l) line emission map, with beam sizes given in Table 2. 



Line 


Ico 


Leo 




Jy km"' 


K km-' pc^ 


'-CO(l-O) 








PdBI 


223 ± 


9 


6.43x10'* 


PdBI + 30m 


485± 


15 


13.5x10** 


i2CO(2-l) 








PdBI 


360± 


11 


2.49x10'* 


PdBI + 30m 


540± 


17 


3.69x10'* 



Table 3. Line intensities and CO luminosities integrated over 
the entire linewidth of the full Av=370 km s ' and the extent 
of the combined PdBI+30m map after correcting corrected for 
primary beam response (typical errors are of the order of 20%). 



The typical 'butterfly' structure visible in the iso- velocity maps 
in Fig. |5]is the signature of a spatially resolved rotating disk. 
In the central 4"-5" the location of the nuclear peak emission 
in emission also moves from the northern part of the 

channel map at the velocity of -170 km s"' to the southeast at 
the velocity of 220 km s"' with respect to the phase-tracking 
center of the map. The displacement of the nuclear component 
as a function of velocity shows that it is resolved, and it im- 
plies the rotation of the nuclear gas around the very center of 
NGC 6574. In the '^CO(l-O) integrated intensity map shown 
in Fig. [3] (left), the molecular gas distribution reveals a regular 



ring-like structure, with a mean angular diameter of ~50" and 
the peak of the emission in the central 4"-5". The integrated 
map was produced by averaging the channel maps over a total 
velocity width of 415 km s"' (from -180 to 235 km s"'). For the 
'^CO(l-O) emission, the ring consists of some molecular cloud 
complexes that appear to be connected to each other They form 
the two spiral arms of NGC 6574 that stait at a PA of ~ 1 00° and 
a projected radius of r ~1.1 kpc (1.55 kpc deprojected with an 
inclination of 45° - see below), near the end points of a promi- 
nent nuclear bar and extend out to radii of about r ~2A kpc 
towards the north and south, approximately along the major 
axis. A comparison to NIR images (Kotilainen et al. 2000) re- 
veals that the position of the central, resolved CO source co- 
incides to within less than 0.25" with that of the AGN. From 
the '^CO(l-O) map in Fig.[3](left) we estimate a gas bar PA of 
about 115° +5°, taking into account that the line emission on 
the western and eastern sides of the bar appear to be displaced 
by about l"-2" to the south and north, respectively. Our value 
of 1 15° is in reasonably agreement with the value of 105° (to- 
wards the western tip of the bar and derived from a lower res- 
olution map) given by Sakamoto et al. (1999ab). It also appar- 
ently agrees with the NIR maps presented by Kotilainen et al. 
(2000). The NIR stellar bar has a PA of 105°, which means that 
our high-resolution CO maps are able to detect a slight advance 
phase shift between the stellar and gas bars, the latter leading 
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by about 10°. Kotilainen et al. (2000) also interpret the elonga- 
tion of the nuclear component as a possible nuclear bar, with 
an orientation of PA= 150°, however, this feature is not certain 
and could be due to projection effects, as the PA of the major 
axis is 165°. The CO spiral arms, starting at the end of the bar, 
mostly lie at the trailing edges of the nuclear bar This particular 
geometry may determine the feeding budget for the gas in this 
region. The spiral arms and the central source are connected 
by a molecular gas bridge - the small bar. Contrary to '^CO(l- 
0), the i2cO(2-l) integrated emission is detected mainly inside 
the central component within the radius of 4"-6", as shown in 
Fig.[3](right). The maximum of the emission is concentrated at 
the center as well. While our '^CO(2-l) map lacks sensitivity 
outside the central 500 pc of the disk, it does give a sharp im- 
age of the molecular gas distribution in the vicinity of the AGN. 
Also the tips of the nuclear bar in the eastern and western parts 
of the galaxy are indicated. 



respectively. In both cases we used equally spaced velocity 
channels with a velocity range from -170 to 220 km s"' (see 
Fig. |5]) with the central channel centered on the systemic ve- 
locity. The maps show the increase in velocities from north to 
south in both lines. The iso-velocity lines of the central 10" are 
parallel to each other but tilted by about 20° with respect to 
those of the outer disk sections of NGC 6574. This is a clear 
indication of the presence of spiral arms at the location of the 
bends between the iso-velocity contours and the influence of a 
barred potential that expresses itself in an overall tilt of the ve- 
locity field with respect to that of the molecular disk. While the 
'^CO(l-O) line shows the velocity field for the whole galaxy 
(disk -I- nucleus), the '^CO(2-l) line emission only allows us 
to derive the velocity field for the central 4" x4" region. Here 
the overall orientation of the iso-velocity lines is similar to the 
'^CO(l-O) but shown at a higher angular resolution - hence the 
influence of the barred potential is more pronounced. 



3.2. Kinematics of tlie molecular gas disk 

The kinematics of the molecular gas are dominated by its mo- 
tion in a rotating disk. In addition we see the influence of spi- 
ral arms and a bar potential. In Fig. |4]we show the position- 
velocity (p-v) plot along the major axis at PA = 165° and along 
the minor axis at PA - 75°. By comparing the 30m and PdBl 
spectra extracted from the corresponding maps at the same po- 
sitions, we find that there is some excess extended emission 
(as also suggested by Fig.[T]i, especially in the southern part of 
NGC 6574. Therefore both the PdBl and the 30m data indicate 
that the molecular gas emission is asymmetric with respect to 
Vsys, so the CO emission is preferentially redshifted (see Fig.[Tli. 
Our results agree well with the previous finding by Sakamoto 
et al. (1999ab), which were obtained at a lower resolution. 

The position velocity diagram along the major axis shows 
the rotation curve that still has to be corrected for the incli- 
nation of the galactic disk. For both lines of '^CO emission, 
the results are similar; however, differences between them can 
be found in the central 2". The p-v diagram of the '^CO(l-O) 
line resolves the disk and nuclear region of the galaxy. There 
is some evidence of '^CO(2-l) line emission at the very cen- 
tral position, although the diagram clearly shows that the emis- 
sion at the systemic velocity is suppressed with respect to the 
emission arising at 0.5" to a 1" separation from the nucleus. 
This expresses itself in two emission clumps neighboring the 
central position and systemic velocity in the '^CO(l-O) p-v di- 
agram. This could hint at the existence of a small nuclear ring. 
The minor axis feature is a clear indication of non circular mo- 
tion, expected in a bar, supporting the nuclear bar hypothesis. 
In the '^CO(2-l) p-v diagram this contrast between the clumps 
is also visible on the major axis and reduced on the minor axis. 
However, the '^CO(2-l) map has less sensitivity, due to its in- 
creased spatial resolution. 

The distribution of the '^CO emission, as well as the iso- 
velocity maps, also indicates a more complex molecular gas 
kinematics in the center of the galaxy. The iso-velocity map of 
'^CO(l-O) and '^CO(2-l) line emission were made from maps 
at flie full angular resolution of 2.48"xl.62" and 0.93"x0.5", 



3.3. Molecular gas masses 

The velocity-integrated '^CO(l-O) flux within the ~50" pri- 
mary beam of the combined IRAM PdBl/30m data is used to 
estimate molecular hydrogen masses. We use /^ot^ ^ ~ 
0.12xlO-^-/co[Jy km s-^]-v-^[Hz]-0-Yl with the observing 
frequency v and the beam size of the corresponding map 6. The 
column density of molecular hydrogen can be derived from /^^ 
by assuming the validity of an average conversion factor X - 
N{H2)irco = 2.2x102" [cm-^K km s"'] (see e.g. Solomon & 
Barrett 1991): 

-2 

cm 

N{H2)^X[— ^]-/^o[Kkms-i] . (1) 

Kkm s ' 

To estimate the molecular gas mass, the column density 
NiHi) then needs to be multiplied with an effective projected 
area A over which the CO line emission is distributed: 

MhAMq] = 1.51 X lO"''* ■ A^(//2)[cm-2] . A[kpc2] . (2) 

Including the relative contribution of helium and and neutral 
atomic hydrogen, the total gas mass contained in the disk of 
NGC 6574 can be estimated via 

Mg,AM^]^\36-{MH,+MH,), (3) 

with helium mass contained in the factor of 1.36. 

The total dynamical mass of NGC 6574 is derived from the 
virial theorem by using the inclination-corrected circular ve- 
locity for a given radius via M^,.„[M0] = 232 x v^Q,(r)[km s"'] 
X r[pc]. The error of the Keplerian dynamical mass due to the 
flatness of the mass distribution is at most 30% for an exponen- 
tial disk (Binney & Tremaine 1987). All results are summarized 
in Table m separated into two components: the nucleus, which 
contains the flux of the inner 4", and the disk, containing the 
residual flux within the 50" of the integrated map. This central 
source contains a molecular gas mass of ~ 10^ Mq. Estimated 
from the kinematics revealed by the '2CO(1-0) line emission, 
the total dynamical mass of NGC 6574 contained within a ra- 
dius of about 15" is M^„,[Mq] ~ 3.5 x 10'°. The molecular gas 
contributes about 6% to the dynamical mass in the inner 6 kpc. 
This is a typical value for a moderately gas-rich barred spiral. 
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Fig. 4. Position-velocity diagrams of '^CO(l-O) and '^CO(2-l) along the major axis at 165° (upper panels) and minor axis at 
75° (lower panels). For the top and bottom diagrams, the left side corresponds to the approximate south and east directions, 
respectively. The color bar gives the relative peak intensities in percent, with 0, 50, and 100% levels indicated just below the 
corresponding marks or the top of the bar. From the top left to the bottom right graph the contour lines are equidistantly spaced 
in 7, 17, 10, and 17% of the peak intensity, respectively. 



3.4. CO line ratios 

Comparison of the two CO line maps, at the same resolution 
and with the same spatial frequency sampling, can yield infor- 
mation about the excitation condition of the gas. We studied the 
variation in the '^CO(2-1)/'^CO(1-0) ratio in areas with signif- 
icant emission levels in both lines, i.e. in channels with veloci- 



ties from -180 to 220 km s"' . Within the 20.5 " primary beam, 
the '^CO(2-l) map was first smoothed to the spatial resolution 
of the '^CO(l-O) map to assure that the two lines sample iden- 
tical regions. The maps were then corrected for primary beam 
attenuation. We find a nuclear '^CO(2-1)/'^CO(1-0) fine ratio 
close to unity and an off-nuclear ratio around 0.2. 
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Fig. 5. Iso-velocity diagram of '^CO(l-O) (left panel) and '^CO(2-l) emission (right panel). The '^CO(2-l) emission is observed 
only in the central part (compare with left panel). The innermost solid contour line is at the systemic velocity and the velocity 
steps are 30 km s"'. 



Parameters 


PdBI 


PdBI+30m 


nucleus: 






IcotolJykms''] 


17.57 


20.30 




1.45 


1.86 


N{H2)[cm-^] 


3.18x10'" 


4.09x10-" 




1.05x10** 


1.18x10* 




1.43x10'* 


1.60x10** 


disk: 








223 


485 


r^gJKkms-^] 


0.83 


2.04 


N(H2)[cm-^] 


1.86x10^" 


4.47x10-" 




6.07x10'' 


1.29x10* 




8.26x10^ 


1.74x10* 




2.04x10'' 


2.13x10' 


dynamical mass: 










3.5 X 10'" 


MH,[Me] 


1.44x10" 


(Mirabel & Sanders, 1988) 



Table 4. Results of IRAM PdBI and PdBI+30m telescope. 



The values are also apparent from the maps in Fig. 3 (see 
caption). The ^^CO(2-1)/'2CO(1-0) ratio is best determined on 



the nucleus with a value close to unity. The deconvolved source 
size in the J=2-l line emission is 1.4+0.3 arcsec. This results in 
a '^CO(2-l) line flux of 55.7 Jy km/s in the 3 mm beam and a 
>2CO(2-1)/12CO(1-0) Une ratio of 0.9+0.1. This ratio is typical 
of moderately dense and warm (n/fj>10'*cm"-', Tj;,„~30 K), op- 
tically thick molecular clouds if a small contribution from the 
cosmic microwave background is taken into account as well 
(see Eckart et al. 1990, Garci'a-Burillo et al. 1993). For the 
weak and extended line emission along the bar, we find an off'- 
nuclear 12CO(2-1)/12CO(1-0) ratio of 0.2±0.1. This is indica- 
tive of a dominant contribution of cold (T^^,,, <10 K) or sub- 
thermally (uht <10"^cm"^) excited molecular gas. 

3.5. Comparison with results at other wavelengths 

It is interesting to compare our results achieved through ob- 
servations of '^CO(l-O) and '^CO(2-l) emission with observa- 
tional results at other wavelengths. Laine et al. (2006) carried 
out 20 cm and 3.5 cm wavelength radio continuum observa- 
tions with the VLA and compare the results with the NIR data 
of H2 and Bry from Kotilainen et al. (2000). In the radio and 
NIR they find a central point-like component and a circum- 
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nuclear ring. Both source components are very consistent with 
our observations. The "jet", which is represented as a connec- 
tion between the nucleus and the circumnuclear ring at 3.5 cm, 
along the major axis of the galaxy, has no counterpart in our 
observations. 

4. Modeling the CO distribution 

In the following section we describe the model used to analyze 
the data obtained with the IRAM telescopes, then derive the 
rotation curve and the possible shape of the gaseous orbits. 

4.1. The model 

The inner 200 pc of NGC 6574 show clear deviations from pure 
circular motions. To analyze the complex kinematics, we mod- 
eled the data with 3DRings similar to the approach of Krips 
et al. (2005), Pott et al. (2004), and Schinnerer et al. (2000b; 
see Appendix B therein). With 3DRings only fully symmetric 
structures can be modeled (i.e., no lopsidedness). The model 
subdivides the disk into many individual (circular or elliptical) 
orbits of molecular gas, which lie on the (possibly tilted) plane. 
With the 3DRings program we investigate the kinematical im- 
print of (warped) rotating gas. 

We model non circular motions with 3DRings via ellipti- 
cal orbits with changing position angles characterizing gas mo- 
tions in a bar The inclination, position angle, and shape of the 
rotation curve for the overall galaxy were held fixed. Each fit- 
ting process was started at large radii and successively extended 
toward the center In each case we tried several start setups that 
all converged to similar (best) solutions with mean deviations 
from the data of less than about 10 km s"' and 0.1 arcsec for 
each radius and velocity in the p-v diagrams and 10° in the 
position angle of the mapped structures. The best decomposi- 
tions were achieved using a systemic velocity of Vsys - 43.75 
km s~' in agreement with the derived position velocity dia- 
grams shown in Fig. |4] (see also spectra and channel maps in 
Figsdl and |2]i. The revised best-fit parameters were obtained 
by comparing the result with the observed PdBI-i-30m data and 
tuning the parameters until the match is optimal. We derived 
an incHnation of 45° +10° and position angle of 165° +10° for 
the galactic disc in space. In Fig.|2]we show the intensity distri- 
bution and rotation curve and in Fig.[8]the ellipticity e and the 
position angle of the elliptical bar orbits as a function of radius 
that we used to model the '^CO(l-O) line emission data. The 
position angle is counted from north to east with respect to the 
northern part of the inclination axis. This is represented by the 
dashed line in Fig.|6] 

As a result of the modeling we find: 

1) The rotation curve is basically flat until a radius of at 
least 1" (this is the resolution limit of the observation and of 
the 3DRings model.) This results in an enclosed mass at this 
radius of about 2.3x10'' Mq. 

2) Starting with almost circular orbits at 20" to 30" radius, 
the ellipticity drops to values of e -0.2 ranging in position an- 
gles from -50° to -90° approaching the bar position angle in 
our data (see also the value of -105° given by Sakamoto et al. 
1999ab). 



3) The model requires a constant disk contribution over the 
entire modeled area with a radius of up to 30". An increased 
central flux component starts at a radius of about 10" and then 
sharply peaks at the center with a value about a factor of 3.0 
higher than the disk component. The FWHM of the nuclear 
component is -1.5"; i.e. it is just resolved, as shown also in the 
maps in Fig.|3] 



200 



-100 




1 2 
Radius in kpc 



Fig. 7. Intensity distribution and rotation curve as a function of 
radius as used for modeling the '-CO(l-O) line emission data. 




1 2 
Radius in kpc 

Fig. 8. The ellipticity e (dashed line) and the position angle 
(full line) of the elliptical bar orbits as a function of radius as 
used for modeling the '^CO(l-O) line emission data. 



5. Computation of the torques 

5.1. Near infrared images 

We used the near-infrared images kindly made available by J. 
Reunanen. The JHK images have been obtained by Kotilainen 
et al. (2000) at UKIRT under 0.6-0.7" seeing and a field of view 
of 70". The H band image is superposed on the CO contours in 
Fig. |9] There is a very good correspondence between the NIR 
image and the CO(l-O) contours, indicating that the molecular 
gas is well aligned along the bar and nuclear bar in NGC 6574. 
However, a slight phase shift can be noticed, with the gas on 
the leading side. The nuclear bar is about 8" (or 1.2kpc) in 
diameter, while the main bar is 20" (or 3.2kpc) in diameter 

The color J-H image is also compared with the CO emis- 
sion in Fig. [To] It reveals a spiral structure emerging of the main 
bar Also, the correlation with the CO distribution is evident. 
The spiral arms in the CO emission wind up in a pseudo-ring, 
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5 -5 
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Fig. 6. Integrated intensity map and position velocity diagram along the major axis for the modeled '^CO(l-O) emission. The 
dashed line at PA = 165° indicates the position of the major axis of NGC 6574. The position-velocity diagram was obtained as a 
cut made through the data cube at this position angle. 



conspicuous in Ha (Gonzalez-Delgado et al. 1997) and in Bry 
(Kotilainen et al. 2000). This star-formation ring is likely to 
correspond to the UHR resonance of the primary bar. Indeed, in 
most studied galaxies, as shown by simulations and confronta- 
tion with observations, the bar ends up near its UHR resonance. 



.a 
I 

I 





Fig. 10. CO(l-O) linear contours (1.7 to 17 by 1.7 Jy/beam) 
superposed on a color image in the near-infrared J-H, from 
Kotilainen et al. (2000), in logarithmic levels. The CO emission 
follows a spiral structure, which winds up in a pseudo ring. 



IC o --0 

Fig. 9. CO(l-O) linear contours (1.7 to 17 by 1.7 Jy/beam) su- 
perposed on the near-infrared H image, from Kotilainen et al. 
(2000), in logarithmic levels. The CO emission is aligned along 
the bar of PA - 105° and also possibly with the nuclear bar, in- 
side a radius of 3". 



5.2. Evaluation of the gravitational potential 

As described in previous papers (e.g. Garcia-Burillo et al. 
2005), we assume that the NIR images can give us the best ap- 
proximation for the total stellar mass distribution because it is 
less affected by dust extinction or by stellar population biases. 
We here recall the essential definitions and assumptions. 

The H image was first deprojected according to the angles 
PA - 165° and / - 45°. We did not deproject the bulge sepa- 
rately, since we do not know its actual flattening, and the galaxy 
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- a late type (Sbc) - does not make a large contribution. The 
image is, however, completed in the vertical dimension by as- 
suming an isothermal plane model with a constant scale height, 
equal to ~l/12th of the radial scale length of the image. The po- 
tential is then derived by a Fourier transform method, assuming 
a constant mass-to-light (M/L) ratio. The M/L value is selected 
to retrieve the observed CO rotation curve. The axisymmetric 
part of the model, fitted by parametric functions, is then derived 
to find the proper frequencies, as shown in Fig.[TT] The rotation 
curve agrees within the uncertainties of the data (Figs. HI with 
the rotation curve shown in Fig. |7] which was derived while 
modeling the CO data alone. 
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Fig. 11. 
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Rotation curve and derived frequencies fi, O - 
for NGC 6574. 
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Fig. 12. Strengths (Qi, Q2, and total Qj-) and phases {<p\ and 
^2) of the m - I and m - 2 Fourier components of the stellar 
potential inside a field of 38" in diameter (r < 3kpc). Note the 
constant phase of the bar. 



For the non-axisymmetric part, the potential <t>(R, 0) is then 
decomposed in the different m-modes: 

^(R, 0) = (i>o(R) + y <:>JR) cos(m6 - <pjR)) 



where ^m(R) and 4>,„{R) represent the amplitude and phase of 
the m-mode. 

The strength of the m-Fourier component, Q,„(R), is defined 
as Qm{R) - ni<i>,„IR\FQ(R)\, i.e. by the ratio between tangen- 
tial and radial forces (e.g. Combes and Sanders, 119811 1. The 
strength of the total non-axisymmetric perturbation is defined 
by 

F'^"\R) 
Qt{R) = — 

where Fj'"'{R) represents the maximum amplitude of the tan- 
gential force and Fo{R) is the mean axisymmetric radial force. 
Fig. [12] displays these values as a function of radius for NGC 
6574. A main bar can be seen clearly, together with an asym- 
metry towards the center. The lopsidedness has a correspon- 
dence in the gas morphology described in Section 3. The bar is 
regular in phase, but its strength is relatively modest. 



5.3. Evaluation of gravity torques 

After having calculated the forces per unit mass (F^ and Fy) 
from the derivatives of <^{R, 0) on each pixel, the torques per 



unit mass t(x, y) can be computed by 

t{x,y) ^ X Fy -y F;, . 

The torque map is oriented according to the sense of rotation 
in the plane of the galaxy. The combination of the torque map 
and the gas density 2(x,y) map then allows derivation of the net 
effect on the gas, at each radius. The gravitational torque map 
weighted by the gas surface density t(x, y) x 2(x, y), normalized 
to its maximum value, is shown in Fig. [13] 

To estimate the radial gas flow induced by the torques, we 
first computed the torque per unit mass averaged over the az- 
imuth, using £(jc,y) as the actual weighting function,i.e. 



t(R)^ 



Jgnx,y)x(xFy-yF,) 



By definition, t(R) represents the time derivative of the spe- 
cific angular momentum L of the gas averaged azimuthally, i.e., 
t(R)-dLldt Ifl. To derive non dimensional quantities, we nor- 
malized this variation in angular momentum per unit time to 
the angular momentum at this radius and to the rotation period. 
We then estimated the efficiency of the gas flow with the aver- 
age fraction of the gas specific angular momentum transferred 
in one rotation (T„,r) by the stellar potential, as a function of 
radius, i.e., by the function AL/L defined as 



AL 
T 



dL 
dt 



1 

X — 
L 



t(R) 



12 



E. Lindt-Kiieg et al.: Moleculai- gas in NUclei of GAlaxies (NUGA) VIII: NGC 6574 



where Lg is assumed to be well represented by its axisymmetric 
estimate, i.e., Lg - R x v,.o,. The AL/L curves for NGC6574, 
derived from the CO(l-O) or the CO(2-l) gas, are displayed in 
Fig.m 




relative R.A. in ai'csec 

Fig. 13. The CO(l-O) contours are overlaid onto the gravita- 
tional torque map (t(x,y) x 2(x,y), as defined in the text) in the 
center of NGC 6574. The deprojected torque map (grey/color 
scale) is normalized to the maximum absolute value of the 
torques. The derived torques change sign as expected in a but- 
terfly diagram, delineating four quadrants. The orientation of 
quadrants follow the bar orientation in NGC 6574. In this de- 
projected picture, the major axis of the galaxy is oriented par- 
allel to the abscissa. The inclined line reproduces the mean ori- 
entation of the bar (PA= 105° on the projected image). 



5.4. Results and discussion 

We show in Fig. [T3]that the derived torques change sign fol- 
lowing a characteristic 2D butterfly pattern. The CO contours 
clearly reveal that the majority of the gas at i~ 0.5-1.5 kpc 
is offset from the bar on the leading edge where the torques 
are negative. The rotation sense in the galaxy is counterclock- 
wise, and the spiral structure is trailing. Towards the center (r< 
400pc), however, the dominating torques are positive. Indeed, 
the average value of dL/L per rotation over 0-400pc is -0.005, 
and over 400- 1 500pc is 0. 1 . 

The interpretation of these results can be done in the clas- 
sical scenario of angular momentum transfer in barred galax- 
ies. The main bar, 1 .6 kpc in radius, ends in a pseudo-ring at 
its ultra-harmonic resonance (UHR), where its pattern speed 
equals Q - k/4. According to Fig. [TT] we can then derive the 
probable pattern speed of ~ 60 km/s/kpc. Its co-rotation should 
then be around 2.8kpc. Since the rotation curve gradient is quite 
strong in the center, as revealed by the CO kinematics, there 
must be an ILR in the circumnuclear region, which might cor- 
respond to the high gas concentration there, at r< 600 pc. The 
ILR is usually the site of high gas concentration in a ring. A 
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Fig. 14. The torque, or more precisely the fraction of the an- 
gular momentum transferred from/to the gas in one rotation 
dL/L, plotted for CO(l-O) (left), and CO(2-l) (right). The 
curve is noisier for the CO(2-l) line due to a lower sensitiv- 
ity at 230 GHz. Both curves, however, show negative torques 
around 1 kpc radius (more precisely from 0.4 to 1.5kpc for the 
CO(l-O)). 

suggestion of such a nuclear ring has been given in Section 
3, by the two prominent concentrations of CO in the major 
axis p-v diagrams. It is obvious that the CO emission is not 
peaked in the very center, around systemic velocity, but corre- 
sponds to an elongated nuclear ring, with highly non-circular 
motions, at a radius of about 130pc. With the pattern speed of 
~ 60km/s/kpc, two ILR must exist, as suggested by the rotation 
curve of Fig. 11. However, this indication is only derived from 
the Q. - k/2 curve with the epicyclic approximation, which is 
not valid for a strong bar. The effective curve is lower, and the 
ILR might be expected between 30 and 900 pc. Moreover, if a 
nuclear ring exists, it is likely to be elongated, corresponding 
to the range of radii usually found in numerical simulations. In 
addition, the near infrared images give a hint of a nuclear bar, 
but this is not certain. This nuclear bar is expected from theory 
to decouple in such conditions, inside the nuclear ILR ring. 

The gas flow towards the center is taken just in the act, at 
the present epoch. The flow rate is not very high, due to a bar 
that is relatively weak (see Fig. [T2] |. But the gas flow already 
has a visible consequence: a high concentration of molecular 
gas towards the center, which must be recent, since the center 
is not the location of any starburst (no peak in Bry or Ha). 

The maps presented here lack spatial resolution to deter- 
mine unambiguously whether the gas is stalled at the ILR of 
the main bar, or if it could be still driven inwards due to a nu- 
clear bar In either case, we conclude that we are seeing the gas 
inflow fueling the center, which soon will feed the AGN nu- 
cleus. A possibility is that the feeding is intermittent, since the 
gas flows induced by gravity torques in the very center (r < 400 
pc) have a very short time scale of tjy„ ~ 12 Myr. 

6. Discussion and summary 

In this paper we present the analysis of subarcsec -resolution 
PdBI observations of molecular gas in the Seyfert 2 NGC 6574 
with a short spacing correction derived from IRAM 30m tele- 
scope observations. Our data complements the analysis made 
by Sakamoto et al. (1999ab) for '^CO(l-O) observations and 
expands the analysis with observations of the '^CO(2-l) line 
emission. The successful descriptive modeling of the disk dy- 
namics with elliptical bar orbits corroborates the hypothesis 
that the molecular gas is strongly influenced by the stellar bar- 
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like potential located within a more extended molecular gas 
disk (Laine et al. 2006, Kotilainen et al. 2000). 

Within the NUGA sample, a whide variety of morpholo- 
gies and dynamics are observed (see introduction and e.g. 
Garcia-Burillo et al. 2005). With the exception of a nuclear 
warp, NGC 6574 shows all of these attributes. But even a nu- 
clear warp cannot be excluded as it only becomes apparent at 
the highest linear resolution (sub- 100 pc), which can only be 
achieved for the closest nuclei (e.g. NGC 1068 and NGC 3227; 
e.g. Schinnerer at al. 2000a). In most cases, there is very lit- 
tle evidence that any of the observed features can be uniquely 
linked to the activity of the nucleus. In some cases (NGC4826: 
Garcia-Burillo et al. 2003, NGC7217: Combes et al. 2004) the 
detailed analysis of the observed source properties even appear 
to rather inhibit than to support fueling of the nuclear region. 

Contrary to the majority of NUGA galaxies studied so far, 
there is clear evidence in NGC6574 of gas inflow on a very 
small scale, of a few hundred parsec, where we find CO emis- 
sion down to our spatial resolution. The gas appears to accumu- 
late at this central location, and the inflow must be recent, since 
there is still no nuclear starburst observed. The CO position- 
velocity diagrams reveal a hint of an elongated nuclear ring 
there, which would be the ILR of the main bar. Our resolution 
does not allow us to go further to determine whether the gas 
is stalled at this ILR ring or driven inwards to feed the AGN 
through a nuclear bar. 

In either case, our result is consistent with the previous 
finding that low-luminosity AGN (the NUGA sample includes 
Seyferts and LINERs) do not require any efficient fueling to 
sustain their luminosity (see e.g. Ho 2003). It is also unclear 
which processes dominate the infaU of gas into the region with 
radii that are beyond our current angular resolution i.e. lin- 
ear scales of less than ~100 pc. Nuclear bars (Shlosman et 
al. 1989), lopsidedness or m = 1 instabilities (Garcia-BuriUo 
et al. 2000), and nuclear spiral density waves (Englmaier & 
Shlosman 2000) provide processes that enable the gas infall 
into the central few parsecs. In the case of NGC 6574, both 
the i^CO(l-O) and the ^^CO(2- 1 ) integrated Une-intensity maps 
consistently show evidence of a weak extension toward the 
southeast. This may be interpreted as indicating lopsidedness 
and may therefore be responsible for transporting molecular 
gas into the central 100 pc. If, however, none of the above- 
mentioned mechanisms are efficient, then as an alternative, vis- 
cosity may support such a transport. 

How viscosity may be responsible for transporting gas to- 
ward the center of such disks has been shown by Duschl et 
al. (2000). The authors suggest that viscosity within these 
gaseous disks may provide an eflftcient AGN fueUng mech- 
anism. However, viscosity could also counteract low-gravity 
torques on the gas. Garcia-Burillo et al. (2005) discuss the de- 
tails of viscosity efl'ects versus gravity torques to drive AGN 
fueling. The authors find that such a counteraction may occur if 
it acts on a nuclear ring of high gas-surface density contrast and 
a few 100 pc size. Whether only one of these mechanisms of 
diff^erent efficiency or a combination of them usually accounts 
for the gas infall is still unknown. Garcia-Burillo et al. (2005) 
propose an evolutionary scenario in which gravity torques and 



viscosity act in concert to produce recurrent episodes of activ- 
ity during the typical lifetime of any galaxy. 

Alternatively, the timescales for fuehng on the few 100 pc 
scale and the onset of activity are so difl'erent that both cannot 
be observed simultaneously (see e.g. Combes 2004). In that 
case detailed investigations of nuclear gas fuehng that drives 
the AGN accretion and possibly the star formation process in 
the very center of the nuclear star cluster have to resort to stud- 
ies of the molecular and atomic gas dynamics on size scales of 
less than 10 pc. This could become possible with the high sen- 
sitivity and angular resolution that wiU be provided by ALMA. 
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